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I. INTRODUCTION 

Hadrons are composed of quarks and anti-quarks and are considered to be governed by 
Quantum Chromodynamics, at least in principle. Since QCD describes a strong coupling 
interaction, perturbative calculation of physical quantities of hadrons is not so reliable other 
than the deep inelastic region where the coupling constant becomes weak due to asymptotic 
freedom and hence other methods like lattice gauge theory have been developed to take into 
account nonperturbative effects. However, the situation has dramatically changed when it 
is discovered that the system of heavy hadrons, composed of one heavy quark Q and light 
quarks q or anti-quarks q, can be systematically expanded in 1/rriQ with a heavy quark 
mass, rriQ. The numerator of this expansion in 1/rriQ could be either Aqcd or m q . 

This theory, HQET (Heavy Quark Effective Theory), Jl| is applied to many aspects of 
high energy theories and many kinds of physical quantities of QCD which can be pertur- 
batively calculated in 1/rriQ. Especially those regarding B meson physics, e.g., the lowest 
order form factor ( which is now called Isgur-Wise function ) of the semileptonic weak decay 
process B — ► Dlv and the Kobayashi-Maskawa matrix element V c b, have been calculated by 
many people. However, since applications of HQET to higher order perturbative calcu- 
lations are very restricted, only forms of higher order operators are obtained. Their Wilson 
coefficients are calculable but some of the matrix elements of those operators are obtained 
so that the whole quantity be somehow fitted with the experimental data. || This is be- 
cause most of the calculations based on HQET do not introduce realistic heavy meson wave 
functions and hence there is no way to determine those quantities completely within the 
model. 

In the former papers |4]|| , using the Foldy- Wouthuysen-Tani transformation || we have 
developed a formulation so that the Schrodinger equation for a Qq bound state can be 
expanded in terms of 1/tuq, i.e., the resulting eigenvalues as well as wave functions are 
obtained order by order in 1/ttiq. In this paper, as one of the applications of our formulation 
we will calculate heavy meson spectrum of D and B, and their higher spin states. In order 



to do so, we would like to start from introducing phenomenological dynamics, i.e., assuming 
Coulomb-like vector and confining scalar potentials to Qq bound states (heavy mesons), 
expand a hamiltonian in 1/tuq then perturbatively solve Schrodinger equation in I/ttiq. 
Angular part of the lowest order wave function is exactly solved. After extracting asymptotic 
forms of the lowest order wave function at both r — > and r — > oo and adopting the 
variational method, we numerically obtain radial part of the trial polynomial wave function 
which is expanded in powers of radial variable r. Then fitting the smallest eigenvalues of a 
hamiltonian with masses of D and D* mesons, a strong coupling a s and other parameters 
included in scalar and vector potentials are determined uniquely. Using parameters obtained 
this way, other mass levels are calculated and compared with the experimental data for D/B 
mesons up to the second order of perturbation. The lowest degenerate eigenvalues of the 
Schrodinger equation gives the so-called A parameters for u, d and s light quarks, which is 
defined by 

where Eh is a calculated heavy meson mass and nig a heavy quark mass. || Meson wave 
functions obtained thereby and expanded in 1/tuq may be used to calculate ordinary form 
factors as well as Isgur-Wise functions and its corrections in 1/rriQ for semileptonic weak 
decay processes. 

All the above calculations are calculated up to I/tuq and analyzed order by order in 
l/rriQ to determine parameters as well as to compare with results of Heavy Quark Effective 



Theory, e.g., the parameters, Ai, A2, and A in Sect. |TV|. The final goal of this approach is to 
obtain higher order corrections to Isgur-Wise functions, decay constants of heavy mesons, 
and the Kobayashi-Maskawa matrix element, V^, by using wave functions of heavy mesons 
obtained so that heavy meson spectrum is fitted with the experimental data. 

Below in Sects. [TT] and [TTI| we will first give formulation of this study and next in Sects. [TV] 
and [V] give quantitative and qualitative discussions on the obtained results. 



II. HAMILTONIAN 



The hamiltonian density for our problem is given by 



H = J dx 3 q ]c (x) (a q ■ p q + (3 q m q ) q c (x) + Q\x) (a Q ■ p Q + (3 Q m Q ) Q(x) , (1) 

Hint = / / dx 3 dx' 3 q c (x) Oi q c (x) V t (x - x') Q(x') O t Q(x'), (2) 

where we consider only a scalar confining potential, O s — 1, V s — S(r), and a vector poten- 
tial, O v = 7 M , V v = V(r), with a relative radial variable r, which we think the best choice to 
phenomenologically describe the meson mass levels. J7|,[| The state of Qq is defined by 

|V) = fd 3 x fd 3 y ^(x - y) q c a \x) QJ,{y) |0) , (3) 



where q c (x) is a charge conjugate field of a light quark q and the conjugate state of Qq by 
(ip\ = \ip) with (0| = |0) . From these definitions, we obtain the Schrodinger equation as 

Htfj = (m Q + E)^, (4) 

where the bound state mass, E, is split into two parts, vhq and E (= E — tuq), so that it 
expresses the fact that the heavy quark mass is dominant in the bound state, Qq, and ip is 
nothing but the wave function which appears in the rhs of Eq. (|3|). 

Operating with the FWT transformation and a charge conjugation operator, which are 
defined in the Appendix A, only on a heavy quark sector in this equation at the center of 
the mass system of a bound state, one can modify the Schrodinger equation given by Eq. (|J) 



(H FWT - m Q ) <8> iPfwt = E ip FWT , (5) 

where a notation £g> is introduced to denote that gamma matrices of a light anti-quark is 
multiplied from left while those of a heavy quark from right. The problem of this paper 
is to solve this equation, Eq. (^) in powers of I/tuq. As described first in this section, 
interaction terms are given by a confining scalar potential and a Coulomb vector potential 
with transverse interaction and a total hamiltonian is given by 



H= {d q ■ p q + f3 q m q ) + (a Q ■ p Q + P Q m Q ) + p q (3 Q S 



^\ l - 2\d q -OiQ + {d q -n){dQ-n)}\V, 



(6) 



where scalar and vector potentials are given by 



Sir) 



+ h V(r) = - 



4 a, 
3 r 



and n 



(7) 



and the vector potential is averaged over longitudinal as well as transverse as given in the 
last term of Eq. (Q). The transformed hamiltonian is expanded in l/rriQ as 



Hfwt — m>Q = H-i + H + Hi + H 2 + ■ 



(8) 



where 



H- 1 = -{l + p Q )m Q , 

H Q = d q -p + (3 q m q - (3 q (3 Q S + 1 1 + - [ d q • 3q + (d q - n) (d Q • n) ] J> V, 
Hi = -- P Q p 2 + —fi q d Q -(p + -qj S + ^-7q • qV 



H, 



2m Q 
1 


2mQ 

2m Q Pq 
1 



m Q 



[ a q + (a q ■ n) n } V, 



(9a) 
(9b) 

(9c) 



1 \ i - 1 i 



Ami 



5m; 



Q 



Ami 



■im. 



(p + q) (a Q -p) +p(a Q - (p + q))+iqxpi Q ] ■ [d q + (a q ■ n) n] V, (9d) 



Here Hi stands for the i-th order expanded hamiltonian and since a bound state is at rest, 



P = Pg = -PQ, P =P g = -Pq , q = P ~P, 



(10) 



are defined, where primed quantities are final momenta and the relation of these momenta 
with particles is depicted in Fig. [I]. 



FIGURES 



Fig. 1 




Definition of momenta of Qq 



FIG. 1. Each momentum is defined. 



Details of derivation of equations in this section are given in the Appendix A. 



III. PERTURBATION 



Using the hamiltonian obtained in the last section, we give in this section the Schrodinger 
equation order by order in I/itlq. Details of the derivation in this section are given in the 
Appendix B. First we introduce projection operators: 

l±0a 



A- 



(11) 



which correspond to positive-/negative-energy projection operators for a heavy quark sector 
at the rest frame of a bound state. These are given by (1 ± f)/2 in the moving frame of a 
bound state with v M the four- velocity of a bound state. Then we expand the mass and wave 
function of a bound state in I/thq as 



E = E - m Q = E e + E[ + E\ + . . . , 

^FWT = ^o + ^1 + ¥2 + ■ ■ ■ , 



(12) 
(13) 



where £ stands for a set of quantum numbers that distinguish independent eigenf unctions of 
the lowest order Schrodinger equation, and a subscript i of E[ and ipf stands for the order 
of 1/rriQ. 
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A. -1st order 

The -1st order Schrodinger equation in 1/mg gives 

whose explicit form is solved in the Appendix C and is given by 

^ = *+=(0 *J m (fO). (15) 

Here £ stands for a set of quantum numbers, j, m, and k and 



1 



M r 



**m(0 = ; . , w . J y*m(n). ( 16 ) 

\-iv k {r)(cr-n)J 

where j is a total angular momentum of a meson, m is its component, k is a quantum 

number which takes only values, k = ±j, ±(j + 1) and 7^ 0, and Ukir) and ffc(r) are 
polynomials of a radial variable r. yj m (Q) are functions of angles and spinors of a total 
angular momentum, j = I + s q + Sq. The corresponding operator for the quantum number 
k is given by 

-/3,(E,./+l), (17) 

which satisfies 

-/? g (E g ./+l)(0 *JJrO) = A;(0 *J m (f) ) , (18) 

i.e., 

'-(3 q (p q -£+l),H^-] =0, (19) 

with /Jq being given in the Appendix D, the lowest order non-trivial hamiltonian, 

Note that since charge conjugation operates on the heavy quark sector the A_ projec- 
tion operator appears in Eq. (|T3|), i.e., positive components of Q corresponds to negative 
components of U C Q. 



B. Zero-th order 



The zero-th order equations are given by 



'a q -p + P q (m q + S) + V]®il% = E<iiJ? l 



0' 



(20) 
(21) 



- 2m Q A + ® tpl + ±A_ [ a q ■ a Q + (a q ■ ft) (a Q ■ ft) ] V ® Vo = °- 

Eq. (^) gives the lowest non-trivial Schrodinger equation with a solution given by Eq. ( |T5*D 
and n = f/r . Detailed analysis of this equation is given in Appendix C. A + components of 
wave functions can be expanded in terms of the eigenfunctions, 



47 = (*} m (f0 o; 



(22) 



Expanding A + <S> ip[ in terms of this set of eigenfunctions, one can obtain the solution for 
Eq. (gJJ) as 



with the coefficients, 



M' 



A+®^ = E c i-^'> 



— ($t'\ [a q -a Q + (a, • ft) (a Q ■ ft) } V \^j) . 



4:m Q 



Here the inner product is defined to be 



(*?|0 



=r 1 



rf 3 rtr(*f(0®^,)), 



and the zero-th order wave functions are normalized to be 1, 



(■qrf *£ \ = 5w5 aP for a, (3 = + or - 



(23) 



(24) 



(25) 



(26) 



C. 1st order 



The 1st order equation is given by 



-2m Q A+ ® ^ 2 + H ® ^[ + H x ® % = Efyi + E[^\ 



(27) 



Multiplying projection operators A± from right with the above equation, and expanding ip[ 



in terms of ^f as 



^ = E(W+ft 



one obtains 



E{ = £ c»' (* 



A+# n A_ 



#7, 



*; 



A_#iA_ 



*; 



(28) 



(29) 



which gives the first order perturbation correction to the mass when one calculates matrix 
elements of the rhs among eigenfunctions and 
1 



jk 



Ei - E k Q 



E 4 e - (* 



A+# A 



#7 ) + ( # 



A_H X A^ 



*t 



c\l = 0. 



for fc ^ £ (30) 
(31) 



This completes the solution for ?/>f since A_^f , or cf_ , is obtained in the last subsection. 
Here we have used the normalization for the total wave function, t/t, as 



i) l \ ijj v ) = 8 U >, 



(32) 



where we have neglected color indices in this paper and hence a color factor, N c = 3, in 
the above equation since it does not change the essential arguments. This definition of 
Eq. fl3"2"|) is admitted because here we are not calculating the absolute value of the form 
factors. The appropriate normalization (normally given by 2E with a bound state mass E) 
will be adopted in future papers in which we will calculate some form factors. This way 
of solving Eq. (|27D is unique and we will use this method below to solve similar equations 
appearing subsection Actually this method has been already used to obtain Eqs. (20 , |21~D 
and to solve Eq. ([H]) obtaining the coefficients c" by Eq. (|24|) . 
One obtains A + <g> ip^ as i n the former subsection, 

A + ®^ = E4-^", (33) 

v 

with the coefficients, 

1 



c 2- 



2mQ 



^„ 



((H -E*)A + ®4 + H 1 A + ®iPt 



(34) 



D. 2nd order 



The 2nd order equation is given by 



-2m Q A+ <g> ip* + H <g> -02 + Hi ® i\)\ + H 2 ® ^o = ^0^2 + #1^1 + ^2^ 



As in the above case (1st order), we obtain 



Ei = E 4i' (* 



A + # n A_ 



*«>+<* 



ffiA_®^) + (* 



A_#,A_ 



m 



which gives the second order perturbation corrections to the mass and 

1. J- 



-2+ p £ _ pk 



E 4 e - (* 



-2+ 



l^/i ,,2 



E 



A + i/nA_ 



*;)-^c?t 



-A' ? 



for fc y^ 



A+#i ®^J 



(35) 



(36) 



(37) 
(38) 



This completes the solution for i\)\ since A_?/>f , or c\i_ , is obtained in the last subsection. 
Although we do not need in this paper, one obtains A + <g> ^3 as 



A+®03 = E C 3-^" 



(39) 



with the coefficients, 

i>i 1 



2m Q 



#7, 



((#0 - ^0) A + ® ^2 + (#1 " *i) A+ ® Vi + ^A+ ® $ 



(40) 



IV. NUMERICAL ANALYSIS 



In this section, we give numerical analysis of our analytical calculations obtained in the 
former sections order by order in 1/ijiq. In order to solve Eq. (^0|), we have to numerically 
obtain a radial part of the wave function, \&^" = ( ^| m ), given by 



*Jm(f) = - 



Uk{r) 



-iVkVr) la -n 



yjmity, 



(41) 
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detailed properties of which are described in the Appendix C. As described in the same 
Appendix, the lowest order, non-trivial Schrodinger equation is reduced into Eq. (|C25| ), 

m , + S + v -a r + l U«.H\ KW\ 

d r + k - -m,-S + V)\v b (r)j °\v k (r))' 

This eigenvalue equation is numerically solved by taking into account the asymptotic be- 
haviors at both r — > and r —>■ oo and the forms of Mfc(r) and v k {r) are given by 

Uk{r), v k (r) ~ w k (r) (-) exp I - (m q + b) r - - ( -J J, (43) 

where 



\/--ra 2 (44) 



4a s 



and Wk(r) is a finite series of a polynomial of r 

iV-l , N i 

«fc(r)=£ «?(-) > (45) 

i=o va/ 

which takes different coefficients for Wfc(r) and Vfc(r). 

i) We have fixed the value of a light quark mass, m q , to be 0.01 GeV as listed in Table | 
since only in the vicinity of this value the D and D* masses can be fitted with the experi- 
mental values. We believe that when the mass, m q = m u = m^, is running with momentum 
these are the correct current quark mass since the momentum is given by the order of the 
B meson mass (~ 5 GeV) |13| though we have not used the running mass to solve the 
Schrodinger equation. We have adopted N — 1 = 7 for the highest power of r that gives 
sixteen solutions to Eq. fl42]) , half of which corresponds to negative energies and another 
half to positive ones of q c state. The lowest eigenvalue of the positive energies is assigned 
to the physical state. That is, although we have a node quantum number, n, other than j, 
m, and k, for Eq. (f|2), we take only the n = solution for each value of k and j quantum 
numbers and we do not take into account higher order node solutions in this paper. The 
lowest positive eigenvalue solution for N = 8 gives the one closest to the zero node compared 
with other N. That is, other N gives a rather oscillatory solution. 
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In the case of a Hydrogen atom, for instance, only the Coulomb potential V ~ 1/r 
survives in the above problem and a radial function, Wk(r), becomes a hypergeometric 
function and its finite series of a polynomial gives discrete energy levels. In our case, since the 
potential includes a scalar term we can not analytically solve the above reduced Schrodinger 
equation, Eq. (^). If we force to make the functions, Uk{r) and Vk(r), finite series and relate 
the coefficients of those functions via recursive equations, it leads us to inconsistency among 
coefficients of each term, r l , of a polynomial. We just assume in this paper that Uk{r) and 
ffc(r) are trial finite polynomial functions of r. 

ii) To first determine the parameters, a s , a, and b appearing in the potentials given by 
Eq. ([?D, and m c , we have calculated the chi square defined by 



x ^(M D -E D f- + ( M D ,-E D .)^ (46) 

°D a D* 

where M DjD * and E DD * are the observed and calculated masses of D and D*, respectively 
and <Jd,d* are the experimental errors for each meson mass. As mentioned already m q = 
m d — m u is fixed to be 10 MeV. Masses, Mp^*, are averaged over charges since we have 
not taken into account the electromagnetic interaction. We have determined the values for 
these parameters which give the least value of x 2 , which is given in Table [Til 

iii) The s quark mass, m s , is determined so that the similar equation to Eq. ( f46|) takes 
the minimum value where we substitute D s and D* meson masses instead of D and D*, 
while the b quark mass, m^, is determined by B and B* meson masses. The input values to 
determine these parameters are given in Table |. 



Table I 



iv) There are two types of solutions to optimal values for these parameters, i.e., one set 
for b < which is listed in Table [TI|, and another for b > 0. However, the solution for b > 
gives large difference between calculated values and observed ones for higher order spins and 
also gives negative values for some spectrum even though the lowest lying sates are in good 
agreement with the observed ones. Hence we disregard this set of parameters. 

Tables [Tl.]| ~ [X| give calculated values, M ca \ c , together with the zero-th order masses, M 
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that are degenerate with the same value of k, ratios, pi/Ma and th/Mq, and the observed 
values, M obs . Here the heavy meson mass, E H , is expanded in I/tuq up to the n-th order as 



it 



e h = m + J2p 1 + Y, ^ ( 47 ) 

with Mq = itiq + Eq being the degenerate mass, pi the i-th order correction from positive 
components of a heavy meson wave function, and n-i the i-th order correction from negative 
components. Note also that the exponential factor in the brackets in the first row of each 
table should be multiplied with a value of each column except for those with the explicit 
exponential factor. 



Table III ~ X 



Strictly speaking each state is classified by two quantum numbers, k and j, and also approx- 
imately classified by the upper component of the light anti-quark sector in terms of the usual 
notation, 2S+1 Lj. Studying the functions yt carefully, one finds the upper component of 
^ m (f) corresponds to the following Table |X|, respectively. 



Table XI 



Here J in J p and 2S+1 Lj is the same as a total angular momentum, j, in the Table [XI 



Although the states can be completely classified in terms of two quantum numbers, k and j, 
we would like ordinarily to classify them in terms of 2S+l Lj. However, the states classified 
by J p = 1 + and 2~, are mixtures of two states in terms of 2S+l Lj as given by the Table 
|Xl|. We would approximately regard the state (k, j) = (l, 1) with " 3 Pi", (-2, 1) with nl P x " 
and (2, 2) with " 3 -D2", respectively, whose legitimacy can be supported by calculating the 
coefficient of each state 2S+1 Lj included in the mixture state. We denote them with double 
quotations so that they remind us an approximate representation of the state in terms of 
2S+1 Lj. Using Eq. ( |C6| ) in Appendix C, their relations are given by 

PWA 1/V5 1W,.*)N 
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(49) 



'" 3 ZV>\ 1 / V3 V2\ f\ 3 D 2 ) 

\ nl D 2 "}) ~ v/5 \-v^ v^J \\ 1 D 2 ) 

v) From Tables pi] ~ |X] we see that the perturbative calculation with these parameters 
might not work well for higher k. Namely masses of 1~ and 2~ for k = +2 give some odd 
values. They become even negative for " 3 _D 2 " of D s and B s at the second order as shown in 
Tables JVUI] and ^. Hence we disregard in this paper all calculated masses of 3 /}i(l~) and 
" 3 I?2"(2 _ ) states in any order. In order to remedy this problem, we may need to improve 
the potential form or adopt some other methods. fl4| Hence here only the first order mass 
spectra are depicted without higher spin states ( 3 -Di(l~) and " 3 D 2 "(2~)) in Figs. ||~|. 

Fig.2 
D meson mass spectrum (first order) 
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FIG. 2. The first order plot of D meson masses. Values in the brackets are the observed values. 
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Fig.3 
D„ meson mass spectrum (first order) 




FIG. 3. The first order plot of D s meson masses. Values in the brackets are the observed values. 
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Fig.4 
B meson mass spectrum (first order) 
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FIG. 4. The first order plot of B meson masses. Values in the brackets are the observed values. 
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Fig.5 
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FIG. 5. The first order plot of B s meson masses. Values in the brackets are the observed values. 



One may also notice that the s quark mass listed in Table |Tj is relatively small (~ 90 
MeV) compared with the conventionally used values, ~ 150 GeV, which is regarded as 
the current quark mass. It is interesting to note that these values are also obtained as 
m s (/x = 2 GeV) in the recent lattice QCD calculations. ]n| 

vi) Two states, pseudoscalar (0~) and vector (1~), are degenerate at the zero-th order in 
1/rriQ since the eigenvalue Eq for these states depends on the same quantum number k = —1, 
which are split into two via the heavy quark spin interaction terms, like —V (a q ■ Eg x n) 
in H 1 and all terms in Hi + given by Eq. ( |U 1 0| ) . Similar resolution of the degeneracy 
among the states with the same value of k occurs via the same interaction terms. 

vii) The simple-minded heavy meson mass formula given by 



m c (E D * - E D ) = m b (E B * - E B ) 



(50) 



holds exactly at the first order calculation. This is because the zero-th order mass of two 
states with the same k is degenerate and by definition the first order correction to this mass 



16 



is proportional to 1/toq as given by H\ of Eq. (||]) or by H" of Eq. ( P10| ). To see Eq. (|50] 
as an prediction, replacing E x with the observed values M x we obtain at the first order, 



M B * - M B m c 

1 L — = 0.299, 



(51) 



M D * - M D m b 

which should be compared with the experimental value, 0.326. This discrepancy between 
the calculated and the observed comes from our calculation of B meson mass spectrum listed 
in Table [V] which give B and B* meson masses slightly different values from the observed 
ones. Hold also similar equations to Eq. ( |50"D for higher spin states with the same k quantum 
number because of the same reason given above. 

viii) The so-called A parameter can be calculated using the definition || 



A = lim (E H - m Q ) = M - m Q = E^ 1 



(52) 



where Eh, M , and uiq are calculated heavy meson mass, the lowest degenerate bound state 
mass, and a heavy quark mass, respectively. Difference of M and tjiq is nothing but the 
lowest leading eigenvalue, Eq, with k = — 1 in our model. From Tables |(| and [TV] and 
m c = 1.457 listed in Table |J one obtains at the first order 



A M = M 0D -m c = M 0D * -m c = 0.412 GeV, 
A s = M 0Ds -m c = M 0Dt -m c = 0.529 GeV, 



(53) 
(54) 



and from Tables [II and IV and m c = 1.347 listed in Table |T|, one obtains at the second 
order 

K,d = M 0D -m c = M 0D * -m c = 0.518 GeV, (55) 

A s = M 0Ds -m c = Mod* -m c = 0.629 GeV, (56) 

where M 0D , M 0D * are the calculated lowest order D meson mass defined by Eq. fl52|). 

ix) Parameters which give nonperturbative corrections to inclusive semileptonic B decays 
are defined as HT6UT 7 ! 
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X t = J- (H(v)\h v (iD) 2 h v \H(v)) , (57) 

A 2 = r^ (H(v)\ h v 9 -a^G^h v \H(v)} , (58) 

ZCLff VTiQ I 

where h v is the heavy quark field in the HQET with velocity v. d B = 3, —1 for pseudoscalar 
or vector mesons, respectively. Then the heavy meson mass can be expanded in terms of 
the heavy quark mass, Ai, A2, and A, as 

z? , a -^1 + d H \ 2 , v 

E H = m Q + A h . . . . (59) 

The first order calculation in 1/ttiq makes 2A2 equal to Eq. ( |50|) and the Ai can be calculated 
using the above equation as 

Ai = 2m b (m b + A U - E B ) - 3A 2 , (60) 

A 2 = ^ (E B . - E B ) , (61) 

where E B and E B , are the calculated B meson masses without the second order corrections. 
The results are given by, at the first order, 

\ 1 = -0.378 GeV 2 , A 2 = 0.112 GeV 2 . (62) 

and at the second order, 

Ai = -0.238 GeV 2 , A 2 = 0.0255 GeV 2 . (63) 

Here we notice that although the first term in Eq. (|60| ) is expected to be ~ 0(1) we find it 
to be small from Table [V] at the first order and obtain the approximate relation, 

Ax ~ -3A 2 . (64) 



These values should be compared with those in Ref. ||18|| , which give A u = 0.39 ± 0.11 GeV, 
Ai = -0.19 ± 0.10 GeV 2 , and A 2 ~ 0.12 GeV 2 . 

x) Recently, it has been pointed out that the kinetic energy of heavy quark in- 
side a heavy meson plays an important role in the determination of the ratios /b//d, 



(Mb* — Mb) / (Md* — Md), and \V u b/V c b\, in which use has been made the Gaussian form 
for the heavy meson wave function and has been adopted the so-called Cornell potential, 
the same as ours. [|19j They have derived the relation of these physical quantities in terms of 



the Fermi momentum, pp, introduced in []20| in which pp is related to a heavy quark recoil 
momentum, p , by 

3 



f 



J d 3 pp 2 <P(p) = -p F . (65) 



where the momentum probability distribution function is given by 

m = te) * e """ y (66) 

They calculated the lhs of Eq. ( |65|) to obtain pp by using the Gaussian form of the wave 



function and then derive the relations between physical quantities and this pp. ]T9[ We have 
the radial wave function given by Eq. (|43|) different from a Gaussian one and hence should 
have relations among physical qunatities and our parameters, a, b, a s , and m c , independent 
of pp and hence we may calculate the lhs of Eq. (^) to check if our calculation gives the 
value similar to other calculations. Our value of the lhs of Eq. (^) gives for {p 2 ) of the B 
meson at the first order 

(p 2 J) = 0.560 GeV 2 , (67) 

and the second order gives {p 2 ) = 0.562 GeV 2 , which should be compared with the latest 
values pp = 0.5 — 0.6 GeV calculated in [19|] which correspond to (p 2 ) = 0.375 — 0.540 GeV 2 . 



xi) When one takes an overall look at the calculated masses, the negative component 
contributions, rii/Mo, are relatively large for both scalar states, 0^, at the first as well as 
second order even though they become very small for higher spin states. Positive components 
constantly contribute to all states. When one compares the first order with the second order 
calculations, one can not conclude that the second order is better than the first as a whole 
although higher spin states of D and B are largely improved at the second order. This 
conclusion may be also supported by the comparison of the first and second order calculations 
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of the parameters, Ai, A2, and A with other calculations. [18|] In order to incorporate the 
second order effects properly, one may need to introduce a different potential and/or method 
from ours as mentioned in v) in this section. 

We have used the following algorithms to numerically calculate the heavy meson masses, 
Gauss-Hermite quadrature to evaluate integrals and the tridiagonal QL implicit method to 
determine the eigenvalues and eigenvectors of a finite dimensional real matrix. |2l] 



V. COMMENTS AND DISCUSSIONS 

In this paper, we have calculated heavy meson masses like -D(s), D?-,, B^, and Bfs based 
on the formulation proposed before, [|J which develops the perturbation potential theory in 
terms of inverse power of a heavy quark mass. The first and second order calculations of 
masses are in good agreement with the experimental data except for the higher spin states 
even though the second order calculation does not much improve the first order. The first 
order calculation of the HQET quantities, Ai, A2, and A, are also in good agreement with 
the other calculations. |18|] A new study on the HQET introduced the Fermi momentum, 
Pf, to obtain other physical quantities. [19|] Although we have not had pp in our mind at 



the begining of this study, the obtained value given by Eq. ( |67D is in good agreement with 
what they have obtained. |19 

We have also found a new symmetry already mentioned in the paper [Q and realized by 
the operator, Eq. (|T7|) , 

-/3,(E,./+l), 

which is always present when one considers a centrally symmetric potential model for two 
particles or when one takes a rest frame limit of a general relativistic form of the wave 
function and is related to a light quark spin structure, i.e., y^ m {VL). That is, this is quite a 
general symmetry, not a special feature peculiar to our potential model. 

One can easily see degeneracy among the lowest lying pseudoscalar and vector states as 
follows. Define 
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|P) = f/- 1 (0 *oo), \V,\) = U£ 1 (0 «7a), (68) 

where the inverse of the charge conjugate operator, £/ c _1 , is defined in the Appendix A 
Eqs. ( |A2|) , *&j m is an eigenfunction obtained in the last chapter, The explicit forms of these 
wave functions are given in the Appendix C, Eqs. ( |C28| , |C29|) . and the quantum number k 
can take only ±j, or ± (j + 1). Assigning these states to D mesons, one can have 

\P) = £>±\ , or D°) , \V,X) = \D*). (69) 



Since these states have the same quantum number k = — 1, these have the same masses as 
well as the same wave functions up to the zero-th order calculation in 1/rriQ. That is, the 
degeneracy among these states is simply the result of the special property of the eigenvalue 
equation. Higher order corrections can be obtained by developing perturbation of energy 
and wave function for each state in terms of Aqcd/^q as given by Eqs. (|i~2| , |13D , 



E = E - m Q = E* + E( + E e 2 + . . . , 

i>FWT = l[>l + i>{ + ^2 + ■ ■ • , 

Finally we would like to discuss qualitative features of form factors/ Isgur-Wise functions. 
Let us think about to calculate form factors for semileptonic decay of B meson into D. 
Taking a simple form for the lowest lying wave function both for B and D as 

where a parameter b is determined by a variational principle, ^(^l/ 15 ^H^/ 13 ) = 0. Then form 
factors are given by 

E(q 2 ) ~ exp const. E 2 (q 2 — g^ ax ) , or £(u) ~ exp const. E 2 {u — 1) 

where 

Q 2 = (pb ~ Pd) 2 , oj = v B ■ v D , g^ ax = (m B - m D ) 2 <-> u max = 1, 

with Vq d being four- velocity of B and/or D meson. This means behavior of form factors 
strongly depends on an eigenvalue, E = E — rriQ of the eigenvalue equation, Eq. ([H]), which 
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is often called "inertia" parameter A q when tuq — > oo. This quantity E does not depend 
on any heavy quark properties at the zero-th order. This result also means that the slope 
at the origin of the Isgur-Wise function includes the term proportional to E 2 . The constant 
term (—1/4 like the Bjorken limit p^fl ) for this slope should be given by a kinematical factor 
multiplied with the above expression. 

To conclude, although there have been various relativistic bound state equations proposed 
so far, nobody has yet determined what the most preferable is. We believe that our approach 
presented here must be a promising candidate. 
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APPENDIX A: SCHRODINGER EQUATION 

In order to derive Eq. (Q) or Eq. ([Sp, we need to calculate the expectation value, 

(i/j\(Ho + H int -E)\ilj), (Al) 

by using the equal-time anti-commutation relations among quark fields, 



l M J X =x' 

{Q*(x),Ql(x')\ , =5 aP 5 3 (x-x'). 



Since the wave function, ip a p(x — y) defined in Eq. (|3]), is normalized to be constant, what 
we need to do is to variate Eq. ( |A1| ) in terms of i> a p(x — y) and to set it equal to zero. Then 
we can easily obtain Eq. (^) with the effective hamiltonian given by Eq. (^|). In the course 
of this derivation, it appears to be clear that p only operates on wave functions while <f only 
on potentials. 
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To derive the FWT and charge conjugate transformed Schrodinger equation given by 
Eq. (H), we have used the following definitions: 

Hfwt = U c Ufwt (p'q) H Up WT (p Q ) E/ - " 1 , ipFwr = U c U FW t (pq) i/j, (A2a) 

Ufwt (p) = exp (W (p) 7 Q • pj = cos W + 7q • p sin W, (A2b) 

$= P -, t^W{p) = ^^ 1 E = ^p 2 + ml (A2c) 

Uc = H° Q l 2 Q = -U; 1 . (A2d) 



Note that the argument of the FWT transformation, Ufwt, operating on a hamiltonian from 
left is different from the right-operating one, since an outgoing momentum, pq ', is different 
from an incoming one, pq. However, here in our study we work in a configuration space 
which means momenta are nothing but the derivative operators and when we write them 
differently, for instance as pq and pq , their expressions are reminders of their momentum 
representation. Hence although the arguments of Ufwt and Up WT look different pq and 
Pq are the same derivative operator, — iV. Here the difference between pq and pq is q 
which operates only on potentials and gives nonvanishing results. The FWT transformation 
is introduced so that a heavy quark inside a heavy meson be treated as a non-relativistic 
color source. The charge conjugation operator, U c , is introduced so that gamma matrices of 
a light anti-quark is multiplied from left while those of a heavy quark from right, which is 
expressed by using a notation ®. 

To derive Eqs. (H), we need to first expand -£/fwt in 1/Vtc.q and then take into account 
the following properties of a charge conjugation operator, U c = i 7q7q, to obtain the final 
expressions, 

U- l Y Q U c = -7S T , (A3) 



i.e., 



U e P Q UZ 1 = -l3 Q T , U c a Q U; l = -a T Ql U C E Q [7" 1 = -2%, 

u c i Q u^ = -%, u clQ u; l = ^ Q T , (A4) 

where a superscript T means its transposed matrix. 
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APPENDIX B: DERIVATION OF PERTURBATION 

In this Appendix we will follow the paper M to derive perturbative Schrodinger equation 
for a hamiltonian given by Eq. (Q) for consistency. Following that paper, we will give an 
equation at each order by using the explicit interaction terms given by Eqs. (H). Here we 
quote the same equations given in the former sections for clarification of derivation. The 
fundamental Schrodinger equation is given by 

(H FWT - m Q ) <g> tppwT = E tp FW T, (Bl) 

and expansion of each quantity in 1/mg is given by 

H FW t - rriQ = #-i + H + E 1 + H 2 + ■ ■ ■ , (B2a) 

E = E e + E{ + E £ 2 + ..., (B2b) 

i>FWT = ^ + ^[ + ^2 + ..., (B2c) 

With a help of projection operators defined by 

a ± = ±4&, < B3 > 

we will derive the Schrodinger equation at each order. 

1. -1st order 

From Eqs. ( pi|) and ( P2|) , the -1st order Schrodinger equation in I/ttlq is given by 

-2m Q A + (8) Vo = 0, (B4) 

which means 

4 = A.® 4. (B5) 

Remember that matrices of a heavy quark should be multiplied from right. That is, the 
zero-th order wave function has only a positive component of the heavy quark sector and is 
given by 
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< = n=\ , (B6) 

where fi and ge are 2 by 2 matrices. More explicit form of this wave function is given in the 
Appendix C. 

2. Zero-th order 

The zero-th order equation is given by 

-2m Q A + ® Vi + Ho ® Vo = #oV>o- (B7) 

Multiplying projection operators, A±, from right, respectively, we obtain 

ff A_®^ = ^Vo» (B8) 

-2m Q A + <g> Vf + # A+ ® ^o = 0, (B9) 

whose explicit forms are given by Eqs. ( p0[) and (|21|), where use has been made of 

A + ® Vo = 0. 



Detailed analysis of Eq. (|B^) is given in the Appendix C. When one expands A + components 

of ip{ in terms of the eigenfunctions, 

ffi(r) 0\ 

«T = _ • (Bio) 

like 

A+<^f = 5>fl'\I/7„ (Bll) 



A_# A + *+) , (B12) 



one can solve Eq. ( |B9| ) to obtain coefficients, c\_ , as 

whose explicit form is given by Eq. (|24]). Here in this paper the eigenfunctions, \l/ £ are 
normalized to be 1, 

(y° tfjA = 5 ee >S a0 for a, (3 = + or - . (B13) 
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3. 1st order 



The 1st order equation is given by 



l„i.l i Trt„i.e 



-2m Q A + ®if>l + H ®iF[ + H 1 ®% = E%ij>i + E{%. 
As in the above case, multiplying projection operators from right, we obtain 

# A_ ® ip[ + i?iA_ ® Vo = ^o A - ® ¥1 + E i^i 
-2m Q A + (8^2 + #oA+ ® ^f + #iA+ <8> ^o = ^o A + ® ^i- 



(B14) 



(B15) 
(B16) 



The first equation, Eq. ( [H 1 5| ) , can be solved like in the ordinary perturbation theory of 
quantum mechanics. First expanding ij)\ in terms of VP^ as 



4=EKn+c{^), 



(B17) 



and next taking the inner product of the whole Eq. ( P15| ) with ( \1/ 



*S 41 + £<£'<* 



A+HnA- 



V„ 



*; 



A_#iA_ 



one obtains 



*; 



?£ J k 



E*c\l + E?5. 



(k- 



(B18) 



where we have used the orthogonality condition, Eq. ( |B13| ), and the lowest Schrodinger 
equation, Eq. (|B~8|), to obtain the first term of the lhs and two terms of the rhs of Eq. ( |B18| ). 
When one sets k = £ in Eq. ( |B18| ), one obtains 



E[ = E A 1 - (*J A + #oA 



%,,) + ($ 



A-HiA. 



n), 



(B19) 



which gives the first order perturbation correction to the mass when one calculates matrix 
elements of the rhs among eigenfunctions, \P fc , like in the ordinary perturbation of quantum 
mechanics. When one sets k ^ I in Eq. (BIS ), one obtains 
1 



£k 

? 1+ = F £ _ F k 



tt' /,Tr+ 

k 



£ 4- (* 



A + iJnA- 



#7, 



*; 



A_FiA_ 



*; 



(B20) 



The coefficient c^ cannot be determined by the above equation, which can be derived by 
calculating a normalization of the total wave function up to the first order, 
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y \tf') = <W, (B21) 

giving 



4$ = 0. (B22) 

This completes the solution for ip\ since A_^f, or c" , is obtained in the last chapter, 
Eq. ( |B12j ). The definition of the normalization, Eq. (P21| ), is already mentioned in the main 
text below Eq. ([32] and hence we will not repeat that argument here. This way of solving 
Eq. ( |B15| ) is unique and we will use this method below to solve similar equations appearing 
the subsection [T il D| as well. 

Eq. ( |B16| ) gives a A + component of ip^ as in the case of A + <g> ip[ in the subsection T UB , 
i.e., setting 

A + ®^ = E C "'^", (B23) 



one obtains coefficients, Cgi , as 



1 



2m Q 



*7i 



H - E e ) A + ® ^ + #! A+ ® ^) ) . (B24) 



4. 2nd order 

The 2nd order equation is given by 

-2m Q A+ ® ^ + H ® Vl + #i ® Vi + #2 ® Vo = ^2 + ^i + ^ 2 Vo- (B25) 

As in the above cases, multiplying projection operators from right we obtain 

H A„ <g> ^ + ff x A_ ® ^f + # 2 A_ ® Vo = ^o A - ® ^2 + E l A ~ ® ^i + #2% ( B26 ) 

- 2m Q A+ ® ^3 + #oA+ ® V 2 + #iA+ <8> Vi + ^A + <g> ip e 

= E e A + ® ^2 + ^fA+ ® ^. (B27) 



Again to solve the first equation, Eq. (|B26| ), first expanding ip^ i n terms of tyf as 
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^ = £(4i'*J' + <£'*«- 



(B28) 



and next taking the inner product of the whole Eq. ( |B26| ) with ( \1/ 



one obtains 



£o fc 4j + £4i'(* 



A+H n A_ 



*7/ 



*: 



#iA_ ®^ 



+ fJA_ff 2 A. 



*; 



-^0 C 2+ + -^1 C l+ + -^2 &fe> 



(B29) 



where we have used the orthogonality condition, Eq. (P13|) , and the lowest Schrodinger 
equation, Eq. (|B8j), to obtain the first term of the lhs and three terms of the rhs of Eq. 
( |B29[ ). When one sets k = I in Eq. ( |B29| ), since cr 1+ = one obtains 



** = £ 



cy_ < ^ 



A+H A 



#,/ 



^ 



HtA- ®ift) + (VJ A_# 2 A 



V't 



(B30) 



which gives the second order perturbation corrections to the mass when one calculates matrix 
elements of the rhs among eigenfunctions. When one sets k ^ I in Eq. ( |B29|) , one obtains 



Ik 

2+ " 


1 

rp£ rpk 


£ 41' 


(^ + 


A+i/oA- 




+ (*£|a 


_# 2 A- 


*< + ) 





*7/) + W H X K_®^{ 



(B31) 



The coefficient c|+ can be derived by calculating a normalization of the total wave function 
up to the second order, Eq. ( P21|) , which gives 



Jfcfc 



1 



-2+ " ■ 2 £ ' r 1 



fcf 



+ 



.,/>•£ 



(B32) 



This completes the solution for i\)\ since A-^f) or c 2- > * s obtained in the last chapter, 
Eq. (JB2§. 

Although we do not use, Eq. ( |B27| ) gives a A + component of i/j^ as in the cases of A + ® ip[ 
and A + ® ip^ given in the subsections |111B| and [TL1C], i.e., setting 



A+®^3 = £c3-^'> 

one obtains coefficients, c" , as given by Eq. (|40|) 



(B33) 
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APPENDIX C: ZERO-TH ORDER SOLUTION 



There have been a couple of trials to solve Eq. (p0|) . JlO],[Il],0] In order to solve the lowest 
non-trivial eigenvalue equation given by Eq. (p0|), i.e., 



[d q -p + (3 q {m q + S) + V]®^ = E^l (CI) 

we summarize and refresh the previous results. First we need to introduce the so-called 
vector spherical harmonics which are defined by 

Yjm = ~ n Yj , (C2) 

Y$ = r V17, (C3) 

\jj{j + 1) 

Y™ = -in x Y<% = , ^ = n x VI™, (C4) 

where Y m are spherical polynomials (or surface harmonics) . These vector spherical harmon- 
ics satisfy the orthogonality condition: 

dQYJi\ny ■ Y™,(Q) = 6 3J ,6 mm ,6 AB , 

where dQ = sin 8d8d(p. This is nothing but a set of eigenfunctions for a spin-1 particle. In 
this paper we need their spinor representation and also need to unitary-transform them to 
obtain the functions, Vj m , as 

aim \ I *j \ I Uj m \ I jm \ 

, k, ;„,- -J =f/ U.r« (C5) 

y j m / \ j m / \ a j m I \ j m I 

where 

1 /v/J + T VJ \ 

u -^n(^ ^)- (C6) 

Yj m (A=L, M, E) are eigenfunctions of j 2 and j z , having the eigenvalues, j(j + 1) and m. 
Parities are assigned as (— ) J+1 , (— ) J , (— )- 7+1 for A=L, M, E, respectively, and Y m has a 
parity (— ) J . Here 
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j = £ + s q + s Q , (C7) 

and s q = a q /2 and Sq = (Jq/2 are spin operators of light anti-quark and heavy quark, 
respectively. y k m are 2 x 2-matrix eigenfunctions of three operators, j 2 , j z , and a ■ £ with 
eigenvalues, j(j + 1), m, and —(k + 1), and satisfy 

~tr (JdSl y h ;,J y* m ) = 5 kk '5 Jj! 5 mm , (C8) 

Here the quantum number k can take only values, ±j, or ±(j + 1), and ^ since 1q does 
not exit. The eigenfunction of the zero-th order hamiltonian depends on three quantum 
numbers, j, m, and k as described next. Functions, YJ m ' and y!j m , have the following 
correspondence to those defined in Ref. as: 



and 



— (i+l) im 7 irra 

7+1 im —7 j'rra 

2/jm ^2/+ , Vjm^V- , 



jm i(-)' i m i(+)' 



(C9) 



(CIO) 



The explicit expressions of the first few y^ m are given by 



Voo = ~h=, 2/oo = — 7t= (o 8 • n) , (Cll) 

V47T V47T 



I. _ Z _m -.1 ? 



y^ = ^=a m , y[ m = -^=(a.n)a-\ (C12) 

V47T V47T 



where we have used 



F ° = ^=, Y? = i\j — oo&9, Y± = Ti\ — sm8e ±i,fi . (C13) 

V47T V 47T V 71 " 



In order to solve Eq. flCl"| ) , one can in general assume the form of the solution as 

4 = ^j = (0 ¥} m (f)), (C14) 



where £ stands for all the quantum numbers, j, m, and k and 

F fe (0 

1 v, 



*LW = y?™(n), (ci5) 



G k (f) 
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Since the effective lowest order hamiltonian does not include the heavy quark matrices, one 
can exclude the symbol <S> from Eq. (|CT|). The form of a radial wave function is in general 
given by 

F k (f)\ _ I fi k (r) - f 2k (r) (a ■ n) 

Gk(f)l \9ik(r) -g 2 k(r)(a-n) 
Substituting this into Eq. (|C1|) , multiplying y* m from left and using the orthogonality 
equation for y* m , Eq. (|C8| ), the simultaneous equations for fa- and gn~ are obtained and 
after some calculations the final form of the wave function is determined to be either 



(C16) 



* 



jm 



fik(r) 
-92 k (r) (a -n) 



iJ j m 



' fik(r) y 



jm 



.92k {r) y i 



jm . 



or 



* 



-hk{r){a-n)\ _ ff2k{r)y jm \ ( hkif) 



jm 



y 



]m 



y~ k 

i)] mi 



(C17) 



(C18) 



9ik[r) J \9ik{r)yj m / \-g lk (r)(a-n) 

where Pauli matrices a are all for a light quark. Because of Eq. (|C18|) , we generally define 
the eigenf unction, ^ k m , given by Eq. (|C17| ). Then the reduced Schrodinger equation is given 
by 



i ( d r + -) Pi + -P2 + (m q + S (r)) p 3 + V (r) 

i ry I ly 



^k(r) 



ES* k (r), 



(C19) 



with 



Here defined also are 



^k-(r) = 



7i*(r)" 
,92k{r). 



(C20) 



Ota 



Pi 



a n 



o q 



1 t 

1 , P2 

10/ \i 



(1 

Pi V, Pq = | = P3 12X2, 

Vo -l 

-i\ (1 

P3 



-1 



(C21) 
(C22) 



Finally introducing the unitary matrix, 

1 



R 



-i 



R- 1 



1 (T 



i 



(C23) 
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one can transform eigenvalue equation as well as eigenfunctions into 
IrR 



1\ k 

i[d r + -)p 1 + -p 2 + (m q + S (r)) p 3 + V (r) 

iy I iy 



1 -R- 1 \<S> k (r) = E k $ k (r), (C24) 



or 



'm q + S + V 



d r + 



-m. 



-d r + * 

S + F 



■Ufe(r-) 



E k n 



u k (r) 

Vk(r) 



(C25) 



with 



$k(r) 



'u k {r) 



Vk(r 



rR^ k (r) 



rfik (r) 
-ir g 2 k (r) 



(C26) 



Then the solution to Eq. fl20|) is given by 

1 ( u k ( r ) \ k _ N 1 t u k ( r ) V'jmiS 1 ) 



* 



jm 



vU^) = - 



(C27) 



-iv k {r){a-n)/ • \iv k (r) y jm (ty 

Throughout the above derivation, use has been made of formulae given in the next Appendix. 

In order to see the spin-flavor symmetry in our case, the explicit form of each lowest 

order wave function is given as follows in the case of J p = 0~, 1~. That these states 

are degenerate can be easily seen from the eigenvalue equation where the eigenvalue Eq 

depends only on the quantum number k and these states have the same value k = — 1. The 

pseudo-scalar state (J p = 0~) is given by 



%'q 



and the vector state (J p = 1 

E em (° ^ 



U-\ ir) 

4nr \q -iu_i(r) {f -a) 



(C28) 



is given by 





w-i (r) 



[€■ a 



(C29) 



'Aitr \o -iv^{r) (r '■ a) 
where use has been made of Eqs. (|C11 , |C12| ). These are transformed into each other via the 
unitary rotation 



7T^ -, 

exp( -e-a 



(C30) 



Here one has to remember that e 2 = — 1 and also that we omit the U~ l operation on the 
wave function for simplicity Similar degeneracy can be seen for a pair of states with the 
same value of k. 
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APPENDIX D: MATRIX ELEMENTS 



In this appendix, we evaluate matrix elements of the rhs of Eqs. (p9| , |36|) to obtain 
mass corrections, E^, E^, those of the rhs of Eqs. (|30| , [37| , J38| ) to obtain the corrections 
of A_ components of wave functions, cf+> c", and to evaluate Eqs. (p^ , |3~iD to obtain 



the corrections of A + components of wave functions ty\ up to the second order (i = 1, 2), 

1 — > c 2— ■ 

Summarizing A_ — A_ and/or A_ — A + matrix elements of the hamiltonian, the following 
equations are obtained. The rhs of Eq. Q2~3] ) is given here again as 

1 



Jk 



2m c 



*I 



H n 



*f) = 7 (*fc | [a q ■ d Q + (a, • n) (a Q • n) ] V |^ + ) , (Dl) 



which requires to calculate the zero-th order matrix elements. Here one must notice the 
following relation, 



jk _ kl* 



1 



2m Q 



% 



H' 



*; 



(D2) 



The coefficient cf * is given by, for £ ^ k, 



Jk 



2m Q E4iV i : fe H-(^if 1 - 



Using Eq. ( [L)2| ), the first order energy correction is given by 



*; 



(D3) 



£f = 2m Q E|c»' 



+ (* 



AT 



*i 



(D4) 



Simplifying Eq. (0), one obtains 
1 



Jk _ \^je'f'k. 



2m Q 



E fe -^ c«-2Ecf (* fc "&S 



*«> + <* 



#r 



*; 



(D5) 



„f A- 



•&0 _ -^0 



2m Q E #'<£.* + £ c»' *+ flr - vj>+ + c «' *+ ^ 



r+- 



#7 



+ $jnr«7-^ 



r>£ /fc 



for £ ^ fc, 



^ = 2m Q E #'<£ + E (<# («? AT" *£) + 41' (*? ^ 



#7 



*; 



#; 



(D6) 
*+) . (D7) 



33 



Although it is apparent that E\ is always real from Eq. ([D4]), it is not clear whether Eq. (P7|) 
is real or not. We will rewrite Eq. ( |JJ7| ) so that reality of E\ is manifest as follows. 



**=££<£.'<*<- 



m 



+ £ 



**<) 4? - K £ K- > 2Re E 4'J (n HT 

v 



*; 



Ami 



+4m Q ReE^V-:"(*r H[ 



*; 



#7 



tt£ 



tt; 



fl", - 



*7* 



t /' 



(D8) 



whose expression is apparently real. 

In the above derivation, we have used the projected hamiltonian at each order, H" , 
which is defined by 



A a HiAp = A a H { A/3 = AaHj = H { A^, 



(D9) 



r af3 



where H i is composed of Dirac matrices, a, (3, E, and 7 . [12 



Hn 



m 



HI 



Hr 



«<? • P + Pq ( m q + S) + V, 

- [ a g ■ a Q + (d q ■ ft) (a Q • ft) } V, 

a q -p + P q (m q - S) + V, 



1 



2mQ 

-lv 

r 

1 

1 



p 2 + V [(a q -p) -i (d q - ft) d r ] -V i (d q -n) + - (d q ■ S Q x nj 
i (dq-n) - - (dq ■ Sq X 



1 
SPq (d Q -p) + - (P q S' + V) (a Q ■ ft) 



1 



1 

Sf3 q {d Q -p) + - {/3q S' ~ V) {d Q • ft) 



-Pa [P 



—^l 2V [(dq -p) -i (dq-n) d r ] (a Q ■ p) - 2iV' (a q ■ ft) (a Q ■ p) 



(DlOa) 
(DlOb) 
(DIOc) 
(DIOd) 

(DlOe) 
(DIOf) 

(DlOg) 

(DIOh) 

(DIOi) 



-iV (a Q ■ ft) [(d q -p)~ l (dq-n) d r ] + -V (dq • i) 7 J j + ^^ 



v \ si (dq ■ ft) (d Q ■ p) 
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+ [(a q • a Q ) - 2 (a q ■ ft) (a Q ■ ft)] d r + - (a q ■ l) 7 q} , (DIOj) 

Ht = H 2 + , (DIOk) 

H^ + = ^r{/3 9 {p+\l) S +\AV - 1 ((3 q S> + V) (t Q •£)}. (D101) 

Use has been made of the following formulae for the gamma matrices, 

(3 Q A± = ±A±, a Q A± = A T a Q , S Q A ± = A±S Q , 7 Q A± = A T 7 Q = =F«qA ± , 
/3 q E q A ± = ±A ± E Q , 1 5 q A ± = A t1q . (Dll) 

Matrix elements of interaction terms among eigenfunctions are calculated below. Degen- 
eracy can be resolved by heavy quark spin-dependent terms which includes 6lq and/or Sq 
dependent terms, i.e., the last terms of H[ and H 2 together with contributions from 
negative components of the wave functions coming from Hq + and H{ + . 

The formulae necessary for calculating the matrix elements are given below when oper- 
ators, (a q -ft), (& q -t), (vq-p), (aq -t),&q- ipq x ft), (a q - a Q ), (aq -ft), and (aq-p), 
operate on the function, y k m (fi) or /(r) y k m (Q). The symbol, <E>, is used in the same meaning 
for 4x4 gamma matrices, i.e., Pauli matrices for a light anti-quark are multiplied from left 
while those for a heavy quark from right. 

(a q ■ ft) <g> y) m = -yj£, (D12) 

(& q -$)®y k jm = -{k + l)y k jm , (D13) 



k + l 



[& q -p)® f (r) y k m = % [d r + / (r) y > 



* + l| , . , ,, , 



-i[d r + / (r) (a q ■ ft) ® y k w (D14) 



(°q ■ t) 



fa ■ i) 



'yj!i +1) \ l /j(2j + 3) 2^/JuTT) \ fy^ +1) 



Vjm J 2J + 1 \2yjj(j + l) -(2j-i)(j + i)J V y] m 

(D15) 
'v&\ _n + 2 \ (y)^ 

yjL; V o -(j-i)J \y~i 
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[a q ■ (Jq x n' ' • • ' 



> 



V]m ) 2J + 1 \\fj(3 + l) 'J J \VjL 

(D16) 

,yjH\ -2i /-(j + i) y/jti + i)\ (yj% +1) ' 

{&q ■ &Q x n) <g> 



aim J 2 3 + 1 \ yjj(j + 1) -j J \ y] 



J H +1) \ 1 / 2 J + 3 -4v5(JTl)\ /^ +1) 



i/Jm / 2J + 1 \-^j(j + l) 2j-l J V l/j m 

j + 1 \ / jr' + l 

yj m \ I iJ j m 

yjm/ \ yj m 



! 



(D17) 



'yj™\ -i / i -2 v /j(j + i)\ ^jl 1 



(D18) 



•</, 



l) 2j + l V-2VjO' + 1) -1 / \vii 



f(r)yf +1) 



f(r)y J 



j in 



dr- 3 - -2^/j(j + l)(d r + m\ ff(r)y] 



2J + 1 V -2^/j(j + l) - J) " (ft- + I 7 1 ) / V /(r) 1/J, 

(D19) 

ifiDyjt 1 ' 

f( r )yjn 

d r + ^ -2y/j(j + 1) (dr + *&) \ (f{r)y] 



2j + l \ -2y/jU + 1) (ft - ¥) ~ (ft- - I 7 i ) / V /(r) ^ 









1. A — A matrix elements 



a. Firsi order terms 



To calculate Eqs. (|29|, p0|), one needs the following A_ — A_ matrix elements, 
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^l,\Hr\yt) = (n\^{p 2 + V[(a q -$-t(a q -n)d r }-V'(a q -Z Q xn)}\yj) 



= itr J d 3 r^'J -±-[p* + V [(a q -p)-i(a q -n) d r ] - V (a q ■ a Q x n)} ® *J m , (D20) 

where the sets of quantum numbers are given by £ = (j, m, k) and £' = (j, m, k'). Some 
simplification occurs because V(r) ~ 1/r hence V' = —V/r. Each matrix element of the 
first order interaction terms is given below. 

-tr fd 3 r *J m V 2 ® *J m = -tr f d 3 r ttjj (g, -p) 2 ® #* m 



dr $£ 



73^, yryfc t 



..2 + ^+12 







-d.+ 



2 1 fc(fc-i) 



$* 



tr / tfV tt* ' V [(3, • $ - i (a q ■ ft) d r ] ® ** 



jm 



(D21) 



rfr $1 



-V (2d r - *±1 



$A 



(D22) 



y(29 r + ^j 

where $fc(r) is defined by Eq. ( |C26| ) in the Appendix C. Some non- vanishing matrix elements 
of the last term of Eq. ( [D2Q| ) are given by 



m, 



v (d q ■ s c 



x n 



#; 



/ -2k 



^T// 



'0 1 



2 , +1 fdr& k V'[ |$ fc forfc= :-(j4 l),j, 



2A- 



it T// 



1 



•3fcil/dr$^'| |$ fc forfc = j + l, -j, 



2^/dr*}V 



'0 1 



2i+i 



1 



$ -(i+i), 



_27i(^l) t 



2J+1 



r 



1 



$ 



j+ij 



and their complex conjugates. 



(D23) 
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b. Second order terms 



VJ, 



To calculate Eqs. (|36| , 37), one needs the following A_ — A_ matrix elements, 

1 



Ho 



*; 



2m% 1 " 



Pa [P + 



itr />r*?' ' ' 



■ jm 



2m^ 



1^ S+Iav+I^'-F') (a Q -/)j|*+), 
|-/3 9 (p+igj^+^A^+l^y-n (c? Q -/)}®^ m . (D24) 



Each matrix element of the second order interaction terms is given below. 



itr y d 3 r tfjjft (p + ^j 5 ® *} m = / rfr $1 



'5+ 
-5- 



$A 



where 



(D25) 



5+ = S 



,2 , fc(*±l) 



■« + 



fit--). 



Since V = — 4a s /3r and A- = — 47r<5 3 (r), we need to calculate 



1 
-1 

2 



-tr / d 3 r < m f Ai ® ^ fc TO = - |$ fc (0)| 2 <J fc>fc ,. 



(D26) 



(D27) 



Nonvanishing matrix elements of the last term of Eq. (|D24j) are given by 



#7/ 



I(/3 ? <?'-V') (« Q .f) |*+) 



f/drSj 



(fc+i)(2fc-i) fy _ v!_ 

2fc+l 







f dr & 



(*-!)(£ + £) 

n (fc-l)(2fc+l) /<?' , V 

u 2fc-l V r r 



$ fc for k = - {j + 1) , j, 



$ fc for k = j + 1, -j, 



^P/**](f-^)( 



1 0' 







(D28) 



$_ 



(i+i)' 



ViO'+i) 



0' 



4^ /*•*,(?+ S) I ,]"•- 



and their complex conjugates. 
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2. A — A + matrix elements 



a. Zero-th order terms 



Among the many A_ — A + components, that of the rhs of Eq. (E3|) is the only matrix 
element to be needed in the later calculations at the zero-th order, which is again given here, 



¥7, 



Ho + 



*t) = {^v \-{a g -a Q + (a q -n) (d Q ■ ft) ] V | tf +) 



= -tr J <Pr*fJ -[d g -a Q + (d q - n) (a Q ■ ft) } V <g> #J, 
Non-vanishing matrix elements are given by 



and 



#7 



(a g • a Q ) V 



#i 







f/dr$l fc | 1 V$ fc forfc= -(j + 1), j, 



2fc+l 






(T 



1 



y$_ 



(j+i)> 



-i^iP/**;! 1 ii<i. J+1 



2j+l 







3+i-, 



(D29) 



(D30) 



V,, 



(a q ■ ft) {d Q ■ ft) V 



f 2ikiJdr*U 



*; 



-1 



1 o 



V$ k for k = -{j + 1), j, 



2j+l J u/ ^-J 



and their complex conjugates. 



-1 



1 



V$_ 



(i+i)> 



(D31) 



6. Firs£ order terms 



The first order A_ — A + matrix element is given by 
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#7, 



Hi 



*; 



*7' 



1 



-5/3, (3 g .p) + i( ) 9,S" + y , )(5o-") 



* 



€ /' 



-tr / d 3 rfy k ' — 



™>q 



-Sp q (a Q .p) + -(p q S' + V')(a Q -ri) 



Non- vanishing matrix elements are given by, for k = — (J + 1), or j, 



V„ 



S(3 q (a Q ■ p) 



*; 







,'dr + ^ 



,^/**[,,s 



d r 







-d r + 



fc-1 



$fc, 



and 



*r# 



09, S' + V) (a Q ■ n) 



9: 



hfdr$t 



S' + V 



2fc+l 



-S' + V 



$fc, 



ViO'+i) 



2.7 






5' + V 







$_ 



0'+i)> 



*L- (D32) 



(D33) 



(D34) 



o -s'+y, 

and their complex conjugates. Note that when one takes complex conjugate, derivative 
operators do not operate on S' and V but on a wave function, $^. 



Matrix elements of interaction terms, H l and iJ 2 > between eigenfunctions (\I/ 



J mJ 



obtained above are summarized below. Degeneracy can be resolved by diagonal as well as off- 
diagonal matrix elements of the last terms oiH 1 and H 2 together with contributions from 
negative components of the wave functions coming from Hq + and H^ + as mentioned earlier. 
Here matrix elements computed by eigenfunctions with the same fc's are called "diagonal" 
and others with different fc's are "off-diagonal". Calculating all the matrix elements from 
the hamiltonian given in the Section |T|, total mass matrix is given by, up to the second order 
in l/m,Q, 
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^E_i+t/_i,o 




















\ 





E-i+U-i,i 














^1,2 











Ei+Uifi 


























£i+t/i,i 


V xU 




















V -2,. 


E^2 + U-2,l 


























£-2 + ^-2,2 











v l-l 














E 2 +U 2 ,l 





I o 




















£■2+^2,2 / 



(D35) 



where 



■" £.-/.• £.-/■• it rri 1 ) i rK 2 ) T/J _ T/( 1 ) ■? i T/( 2 ) 3 (D36") 



£ fc = m Q + ^ + 2# + E k 2 , U k , j = U^ + Ul% VI k , = Vgj/ + Vt%? 



Here superscripts mean the order of 1/tjiq, k and k' stand for k quantum number, and j 

for the total angular momentum. For instance, V fc J means the matrix element of the third 

,k t 



term of H 2 between ^f TO and ^f m given by Eq. ( P24|) 



As for the A_ — A + matrix element, we only need non- vanishing matrix elements of H { 



-+ 



up to the second order in 1/rriQ, which are given 



or cjk as one can see from Eqs. ([Dl] 
below. 

/ 


£.-» 






V o 

where the integer in the brackets is a value of a total angular momentum, j, and it turns 
out that this matrix is hermitian, i.e., cf^(jo) = c\t{ja)* for the same value of j = jo- 






cr_ M (o) 














\ 








-1,1/11 


<h- C 1 ) 
































:*(!) 














1,2/iS 

CiL(l) 





-*(i) 














C1 -_ 2 ' 2 (l) 























cr_' 2 (2) 








2,l/i\ 
Ci'_(l) 


c?L" 2 (l) 























c?l~ 2 (2) 





o ) 



(D37) 
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TABLES 





TABLE : 


.. Input values to determine parameters (Units are 


in GeV.) 






m q = m u = m d 


M D 


M D * 


M Da 


M Di 


M B 




M B * 


0.01 


1.867 


2.008 


1.969 


2.112 


5.279 




5.325 


TABLE II 


Most optimal values 


of paramet 


3rs determined by the least 


chi square 


method. 


parameters 


ot s 


a (GeV- 1 ) 


b (GeV) 


m c (GeV) 


m s (GeV) m b (GeV) 




x 2 


first order 


0.3998 


2.140 


-0.04798 


1.457 


0.09472 


4.870 


2.5 


x 10~ 12 


second order 


0.2834 


1.974 


-0.07031 


1.347 


0.08988 


4.753 


6.8 


x 10~ 9 






TABLE III. 


D meson mass spectrun 


i (first order) 









State (J p 


M 


Pi /M 


rai/Mo 


M C3i \ c 


M ohs 


% (o- 


) 1.869 


-2.01 (xl0~ 2 ) 


1.93 (xlO^ 2 ) 


1.867 


1.867 


3 s, (l- 




7.37 


0.101 


2.008 


2.008 


3 Po (0+ 


) 2.276 


-0.373 


1.59 


2.304 


- 


" 3 Pi" (1+ 




7.50 


0.0883 


2.449 


2.422(?) 


nl Pi" (1+ 


) 2.216 


7.51 


0.0210 


2.383 


- 


3 P 2 (2+ 




9.54 


4.72xl0~ 7 


2.428 


2.459 


3 Di (1- 


) 2.440 


181 


0.0242 


6.850 


- 


" 3 J D 2 " (2" 




177 


4.28xl0~ 7 


6.747 


- 
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TABLE IV. D s meson mass spectrum (first order) 



state (J 1 



Mn 



Pi /M 



ui/Mq 



M calc 



M obs 



So 


(o- 


') 


3 5x 


(1- 


") 


3 Po 


(0+) 


))3 p » 


(1 + ) 


nl Pl" 


(1 + ) 


3 P 2 


(2+) 


3 A 


(1- 


") 


" 3 D 2 " 


(2- 


") 



1.986 



2.288 



2.335 



2.401 



2.66 (xl(T 2 ) 


1.67 (xl(T 


- 2 ) 


1.966 


1.969 


6.95 


0.0851 




2.125 


2.112(?) 


0.769 


1.45 




2.339 


- 


8.62 


0.0835 




2.487 


2.535 


6.87 


0.0194 




2.496 


- 


8.79 


8.87x10" 


•7 


2.540 


2.574(?) 


159 


0.0282 




6.230 


- 


126 


8.63x10" 


•7 


5.428 


_ 



TABLE V. B meson mass spectrum (first order) 



state (J 1 



M 



Pi /M 



ni/Mo 



M cak 



M obs 



% (0-) 

3 Si (1-) 

3 Po (0+) 
" 3 Pi" (1+) 
»i Pl » (!+) 

3 P 2 (2+) 



3 Di (1- 



n3 D 2 " (2" 



5.281 



5.689 



5.629 



5.853 



■2.13 (xl0~ 3 ) 


2.05 (xlO" 


" 3 ) 


5.281 


5.279 


7.80 


0.107 




5.323 


5.325 


-0.447 


1.90 




5.697 


- 


8.98 


0.106 




5.740 


- 


8.85 


0.0247 




5.679 


- 


11.2 


5.56x10" 


-8 


5.692 


- 


225 


0.0302 




7.172 


- 


220 


5.34x10" 


-8 


7.141 


. 
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TABLE VI. B s meson mass spectrum (first order) 



state (J 1 



M n 



Pi/Mq 



m/Afo 



M cak 



M obs 



'So (0- 



3 Si (i- 



3 Po (0+) 
" 3 Pi" (1+) 
»i Pl » (!+) 

3 P 2 (2+) 



3 I>i fl" 



" 3 ZV (2" 



5.399 



5.701 



5.748 



5.814 



2.92 (xl0~ 3 ) 


1.84 (xlO" 


" 3 ) 


5.393 


7.65 


0.0936 




5.440 


0.923 


1.74 




5.716 


10.3 


0.100 




5.760 


8.35 


0.0236 




5.796 


10.7 


1.08x10" 


7 


5.809 


197 


0.0348 




6.959 


156 


1.07x10" 


7 


6.719 



5.369 



TABLE VII. D meson mass spectrum (second order) 



state (J p 


) M 


pi /M 


m/Mo 


p 2 /M 


n 2 /M 


M cak 


M obs 


% (o- 


) 1.865 


0.457 (xl0~ 4 ) 


1.03 (xl0~ 2 ) 


-1.86(xl0~ 17 ) 


-1.40 (xl0~ 2 ) 


1.867 


1.867 


3 s 1 (l- 




7.41 


0.0540 


-2.95xl0~ 4 


0.209 


2.008 


2.008 


3 Po (0+ 


) 2.319 


2.33 


0.828 


3.44 


0.689 


2.408 


- 


" 3 Pi" (1+ 




7.00 


0.164 


3.63xl0~ 4 


-0.0107 


2.486 


2.422(1) 


nl Pi" (1+ 


) 2.205 


8.22 


0.0532 


-3.82xl0~ 4 


-0.0541 


2.388 


- 


3 P 2 (2+ 




8.67 


0.0156 


-0.157 


0.112 


2.399 


2.459 


3 Di (1- 


) 2.584 


8.29 


0.154 


2.14xl0~ 4 


-0.177 


2.798 


- 


" 3 L> 2 " (2" 




8.13 


0.0133 


3.17 


-0.134 


2.791 


- 
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TABLE VIII. D s meson mass spectrum (second order) 



State (J p 


M 


Pi /M 


m/Mo 


P2/M 


n 2 /M 


M cak 


M obs 


l s Q (o- 


) 1.976 


0.150 (xl0~ 2 ) 


90.6 (xl0~ 4 ) 


2.81 (xl0~ 17 ) 


-1.60 (xl0~ 2 ) 


1.965 


1.969 


3 Si (1- 




7.76 


4.65 


-7.02xl0~ 4 


0.243 


2.133 


2.112(?) 


3 Po (0+ 


2.331 


3.20 


76.8 


2.08 


0.956 


2.446 


- 


" 3 Pi" (1+ 




8.12 


3.87 


2.20xl0~ 3 


0.0361 


2.527 


2.535 


nl Pi" (1+ 


) 2.317 


7.29 


0.925 


-2.22xl0- 3 


0.0229 


2.482 


- 


3 P 2 (2+ 




8.31 


2.96xl0~ 8 


-5.98 


-0.00166 


2.509 


2.574(?) 


3 A (1- 


) 2.582 


702 


0.770 


5.38xl0~ 4 


-0.00632 


2.072 


- 


" 3 L> 2 " (2" 




-4470 


2.65xl0~ 8 


212 


0.000139 


-113 


- 



TABLE IX. B meson mass spectrum (second order) 



State {J p 


M 


Pi /M 


ni/Mo 


P2/M0 


n 2 /M 


M calc 


M obs 


l s Q (o- 


) 5.271 


0.212 (xl0~ 2 ) 


10.3 (xl0~ 4 ) 


2.80 (xl0~ 17 ) 


-3.97 (xl0~ 4 ) 


5.286 


5.279 


3 5x (1- 




0.862 


0.541 


-8.41 xl0~ 6 


0.593 


5.317 


5.325 


3 Po (0+ 


5.725 


0.387 


9.50 


-3.54 


2.24 


5.754 


- 


" 3 Pi" (1+ 




0.965 


1.89 


1.18xl0~ 5 


-0.0347 


5.782 


- 


»ip l » (j+ 


5.611 


1.08 


0.593 


-1.21xl0~ 5 


0.171 


5.672 


- 


3 P 2 (2+ 




1.23 


0.174 


-1.06 


0.352 


5.680 


- 


3 D 1 (1- 


5.990 


1.18 


1.88 


7.40xl0~ 6 


-0.614 


6.061 


- 


" 3 L> 2 " (2" 




1.28 


0.163 


-6.48 


-0.465 


6.066 


- 
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TABLE X. B s meson mass spectrum (second order) 



State (J p 


M 


Pi /M 


ni/M 


P2/M 


n 2 /M 


M calc 


M ohs 


% (o- 


) 5.382 


0.161 (xl(T 2 ) 


9.43 (xl0~ 4 ) 


-7.46 (xl0~ 17 ) 


-4.72 (xl0~ 4 ) 


5.393 


5.369 


3 Si (1- 




0.882 


0.483 


-2.07xl0~ 5 


0.717 


5.430 


- 


3 P (0+ 


5.739 


0.505 


8.84 


3.12xl0~ 4 


3.12 


5.773 


- 


" 3 Pi" (1+ 




1.11 


0.446 


7.19xl0~ 5 


0.118 


5.801 


- 


nl Pi" (1+ 


) 5.723 


1.03 


0.106 


-7.21xl0~ 5 


0.0743 


5.782 


- 


3 P 2 (2+ 




1.19 


3.39xl0~ 9 


-3.78 


-0.00540 


5.791 


- 


3 A (1- 


) 5.988 


43.9 


0.0940 


1.86xl0~ 5 


-0.0219 


8.618 


- 


" 3 L> 2 " (2" 




-136 


3.24xl0~ 9 


-2.06 


4.82xl0~ 8 


-2.174 


- 



TABLE XI. States classified by various quantum numbers. 



k 


-1 


-1 


1 


1 


-2 


-2 


2 


2 


j 





1 





1 


1 


2 


1 


2 


J p 


o - 


l - 


0+ 


1+ 


1+ 


2+ 


l - 


2 - 


2S + l Lj 


% 


3 ^ 


3 P 


3 Pi, l Pi 


^i, 3 Pi 


3 P2 


3 D, 


3 D 2 , l D 2 
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